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ABSTRACT 

We update a physically-motivated model of radiation damage in the Hubble Space 
Telescope Advanced Camera for Surveys/ Wide Field Channel, using data up to mid 
2010. We find that Charge Transfer Inefficiency increased dramatically before shuttle 
Servicing Mission 4, with ~ 1.3 charge traps now present per pixel. During detec- 
tor readout, charge traps spuriously drag electrons behind all astronomical sources, 
degrading image quality in a way that affects object photometry, astrometry and mor- 
phology. Our detector readout model is robust to changes in operating temperature 
and background level, and can be used to iteratively remove the trailing by pushing 
electrons back to where they belong. The result is data taken in mid-2010 that recovers 
the quality of imaging obtained within the first six months of orbital operations. 
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1 INTRODUCTION 

The harsh radiation environment above the Earth's atmo- 
sphere gradually degrades all electronic equipment, includ- 
ing the sensitive Charge- Coupled Device (CCD) imaging de- 
tectors used in the Hubble Space Telescope (HST) Advanced 
Camera for Surveys (ACS)/Wide Field Channel (WFC). 
The detectors work by collecting photoelectrons in a po- 
tential well at each pixel. At the end of an exposure, these 
electrons are then transferred, row by row, to an amplifier 
at the edge of the device, where they are counted. However, 
radiation damage to the silicon lattice creates charge traps 
that temporarily capture electrons and release them only af- 
ter a characteristic delay. Any electrons captured during the 
transfer to the readout register can reemerge, several pix- 
els later, as a spurious "Charge Transfer Inefficiency" (CTI) 
trail behind every bright source. 

CTI trailing is particularly troublesome because the 
amount of flux trailed is a nonlinear function of the flux, 
size and shape of a source. The effect is therefore not a 
convolution. A multitude of ad-hoc schemes have been in- 
vented to estimate (and subtract) the effect of CTI from 
catalogues of object ph o tometry, astrometry a nd shape {e.g. 
Chiaberge et all 120091 ICawlev et all 120021 : iRhodes et d\ 
20071 ). However, since CTI moves electrons around fairly 
predictably at the image level, the ideal approach for correc- 
tion is to directly shuffle those electrons back to where they 
belong. Since detector readout is the last process to hap- 
pen during data acquisition, this can be conveniently carried 
out as the first process in data processing. While no algo- 
rithm can undo the nonlinear movement of electrons in a sin- 



gle step. iBristowl (|2003al ) and lPiatek et al\ l|2005l ) proposed 
an iterative algorithm to remove trailing by repeatedly (re- 
) adding new trailing. This requires a model of the (forward) 
readout process. A ph ysically-motivated m odel was devel- 
oped for ACS/ WFC bv lMassev efall (|2010h . using measure- 
ments from trails behind w arm pixels in science imaging. 

This letter updates the lMassev et al\ l|2010ri CCD read- 
out model and pixel-based CTI correction. In fJ5] we account 
for an additional species of charge trap with a long charac- 
teristic release time, and measure the density of traps in the 
detector up to mid-2010. In 331 we implement the improved 
CTI correction algorithm and evaluate its performance. In 
21 we discuss the overall performance of the detectors, in 
light of changes to their operational temperature and the 
long period during which they were offline. 



2 UPDATED CCD READOUT MODEL 
2.1 Well filling rate 

We use a "volume- driven" CCD readout model, whose first 
ingredient is the rate at which electrons fill up the poten- 
tial well in a pixel. A cloud of n e electrons grows in size 
as electrons are added, and a cloud with a larger cross- 
sectional area will be exposed to more charge traps when 
it is swept thro ugh the silicon lattice during reado ut. As 
first suggested bv lBiretta fc Kozhurina-Plataisl (|2005l ) , if the 
traps are uniformly distributed in 3D, the well filling rate 
can be measured using hot pixels (which would appear as 
isolated 5-functions in the absence of radiation damage) in 
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ordinary, on-orbit imaging - from the increasing fraction 
of electrons trailed behind increasingly warm pixels. Uni- 
versally, the fraction of trailed electrons is greatest for faint 
sources, demonstrating that the size of the cloud grows more 
slowly than the number of electrons. This crucial point ex- 
plai ns why CTI is nonlin ear. 

IMassev et al\ l|2010l ) parameterised the height h(n e ) of 
a cloud as zero for the first d m 100 electrons, then increas- 
ing as ft oc (n e — d) a where a ~ 0.58. The first electrons 
were assumed to reside in a supplementary buried channel or 
"notch" specifically intended to compress their volume and 
minimise the degradation of very faint sources. The notch 
is created by doping the silicon lattice, but HST engineers 
now believe that the initial atomic implant in the ACS/WFC 
detectors was unstable and has diffused (Linda Smith, priv. 
comm. 2010). If this were true, it would result in trailing 
behind even faint sources. The interpretation of our (up- 
dated) measurements of faint in science images is hindered 
by the zodiaca l sky background, but the data are consistent 
with no notch. [Anderson fc B edin (2010) studied warm pix- 
els in dark exposures, which have less noise, and confirmed 
trails behind warm pixels containing as few as ~20 elec- 
trons. We shall therefore adopt a model in which the notch 
is no longer operational d = 0. Strictly, we should model 
the gradual disappearance of the notch over time - but the 
same sky background that makes it difficult to measure this 
effect also hides real science data from its influence. 

The pr ofiles of CTI trails i n dark exposures are repro- 
duced from lAnderson fe Bedinl (|2010T l in figure [TJ The rela- 
tive fraction of electrons trailed behind increasingly warm 
pixels confirms that a ~ 0.57. However, when ignoring 
a notch, mea s ureme nts from science images like those in 
IMassev et all f|2010l s l (see El prefer a = 0.465 ± 0.016 
both before and after Servicing Mission 4. The effect of 
this parameter is apparent as the difference between the 
data and the solid g rey curves of figure [TJ which show the 
IMassev et al\ (|2010T ) prediction but using the lower value 
of a. For very hot pixels, the fi rst ~10 pixels of the pre - 
dicted trail are within 3% of the lAnderson fc Bedinl (2010) 
data (c.f. the best-fit dashed curve) - an impressive agree- 
ment considering these measurements are completely inde- 
pendent. However, predictions of the relative trail heights 
begin to disagree when the model is extrapolated down to- 
wards faint trails. Since the measurements of these faint 
trails are affected by a complex interaction with the sky 
background, we adopt our measurement of a = 0.465 be- 
cause its origin is closest to the data we will eventually want 
to correct. 

2.2 Charge trap species 

The second ingredient of a CCD readout model is the density 
and characteristic release times of charge tra ps. Shockley- 
Read fc Hall theory of solid-state devices (e.g. iHardvIl 19981 : 
iJanesick et al J200ll ) suggests that we can expect several dis- 
tinct species of traps at a variety of energies AE below the 
band gap, all of which capture charge almost instantly then 
release it with a probabilistic delay governed by an exponen- 
tial e~*/ T . The characteristic release time r depends upon 

operating temperature as r oc j"- 2 e AE / k T m 

In early ACS data, IMassev et al . (2010) found two 
species of traps with characteristic release times r = 




Pixels i away from warm pixel 

Figure 1. CTI trailing behind warm pixels in dark expo- 
su res. Black points r e produ ce the measurements from figure 5 
of lAnderson &; B edin ( 2010J) , for warm pixels at least 1500 pix- 
els from the readout register and containing approximately 20, 
200, 2000 and 20000 el ectrons (bottom to t op). Solid grey lines 
show predictions of the IMassev et al. I d2010h model, which used 
only the first 9 pixels behind warm pixels, and has had to be ex- 
trapolated down to the lower trails. The agreement between these 
completely independent analyses is impressive. Dashed lines show 
best-fit models of the trail behind the brightest curve (which is 
most reliably measured), using a double exponential (with the 
decay times of lMassev et aZ.ll2010h and a triple exponential with 
free parameters. The dotted lines add secondary electron capture 
to the calculation, whereby trailed electrons can be recaptured 
and retrailed. This represents the difTerence be tween our full algo- 
rithm and the (much faster) approximation of Ande rson fc Bedinl 



{0.88, 10.4} multiplied by the 3212/us CCD clocking speed, 
and associated them w ith impure E-cen tre complexes at 
AE = 0.31 and 0.34 eV l|Hopkinsonll200 ll ). The trap species 
were present in a density ratio of 1 : 3. 

In July 2006, the operating temperature of the 
WFC detectors was lower e d from -77C to -81C 
dSirianni. Gilliland fc Sembackl 120061 : iMack et al\ 120071 ). 
lAnderson fc Bedinl (|2010t l modelled the trail profiles in 
subsequent imaging using an empirical look-up table. 
However, the smooth curves overlaid on their data in 
figure [TJ demonstrate that the profiles can still be accurately 
fit using multip le exponentia ls . The steep dashed line shows 
a two-species IMassev et al ] <|2010h prediction. The more 
extended dashed line shows the best-fit three-species model 
in which both the trap densities and release times are 
allowed to vary. This analysis yields trap release times of 
r = {0.74T0.55, 7.7±4.3, 37±33} x 3212^s with amplitudes 
of {0.18 ± 0.10, 0.61 ± 0.3, 0.51 ± 0.26} traps exposed to the 
20,000 electron charge cloud - i.e. a ratio of 1 : 3.38 : 2.85. 

Pure E-centre complexes 0.44 eV below the conduction 
band are expecte d to produc e the next-longest trails, with 
r ~ 180x3212ps (|Jonesll2 000). This is longer than measured, 
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but the discrepancy may be due to degeneracies in the fit- 
ting of decaying exponentials. This is a notoriously difficult 
task because the exponentials become more and more sim- 
ilar to a constant as they get longer. Our measurement of 
long trails is even more difficult because the sky background 
in noisy and its subtraction is uncertain. However, the de- 
generacy of exponentials also means that our trails can be 
successfully fitted with different r (and even an additional 
trap species) by simply adjusting their normalisations. We 
therefore adopt the superior three-trap model, with the fit- 
ted parameters after July 2006 and, assuming that the third 
species are indeed pure E-centres, we update the pre- July 
2006 values to r = {0.48, 4.86, 20.6} x 3212tis. 



2.3 Charge trap density 

As in iMassev et ai\ (^OIOT ). we measure the effective trap 
density from the amount of trailing in the first nine pixels 
behind warm pixels in archival HST imaging. Large extra- 
galactic surveys prove most useful to build up a uniform 
dataset extending over a long period of time, and to isolate 
the warm pixels from crowding by real astronomical sources. 
To span the entire lifetime of ACS, we gather data from HST 
programmes GO-9075 and GO-10496 (PI: Saul Perlmutter), 
GO-9822 and GO- 10092 (PI: Nick Scoville), GO- 10896 (PI: 
Paul Kalas), GO-11563 (PI: Garth Illingworth) and GO- 
11600 (PI: Benjamin Weiner). The first four of these were 
observed using a commanded gain setting of 1 and the last 
three with a setting of 2; these are distinguished as grey and 
black data points respectively in figure[2] The exposure time 
and filters - hence the background level - also vary between 
programmes. 

Up to July 2006, we find that the total effective charge 
trap density increases linearly over tim^E The trap density 
extrapolates back to a value at launch of a remarkably low 
value of p q — 0.014 ± 0.04 per pixel. These manufactur- 
ing and process traps were dominated within less than two 
months by radiation-induced traps, which were created in 
orbit at a rate of (3.60 ± 0.26) x 10" 4 per day. 

After July 2006, the traps themselves continued to ac- 
cumulate at a rate of (4.77 ±2.76) x 10 -4 per day. Lowering 
the operating temperature and lengthening the trails imme- 
diately reduced the amount of spurious flux in the first nine 
pixels by 22%. According to our model, however, the same 
amount of flux was lost, but it was just moved further. The 
continuity of the apparent trap density around this time, 
through changes in operating temperature, background level 
and default gain settings, provide a strong vindication of our 
model. Overall, the total trap density until January 2007 is 
well fit by p q = (0.50±0.018)±(t-1359) (3.55±0.22) x 10" 4 
per pixel, where t is the number of days after 1 March 2002. 

Since January 2007, degradation has been more rapid. 
Because of the long period when ACS was offline, it is not 
clear exactly when this damage accrued. It is quite possible 
that the radiation exposure simply increased. The solar cy- 
cle maximum ended around 2006, and the density of charged 
particles in low Earth orbit counterintuitively in creases dur- 
ing solar minimum l|Sirianni fc Mutchlerl [2006) . However, 



even though our data suggest the rate of trap creation in- 
creased slightly in late 2006, it appears to have slowed again 
since 2009. A more likely scenario is that the damage built 
up abruptly while ACS was offline and warm, before shuttle 
Servicing Mission 4. The subsequent trap density is best-fit 
by Pq = (1.25 ± 0.020) + (t- 2873) (2.93 ±2.25) x 10 _4 /pixel. 



iMassev et al demonstrated that monthly annealings do 

not reduce the effective trap density, so we also ignore them here. 



2.4 Algorithmic development 

lAnderson fc Bedinl (|2010l l and IShort et al\ (|201Ch inge- 
niously invoke a first-order symmetry of the readout process 
to increase the speed of the readout (and correction) al- 
gorithm. Electrons beginning 2048 pixels from the readout 
register undergo this many pixel-to-pixel transfers during 
readout. Each time, electrons may be captured or released 
by charge traps but, if the number of free electrons is high 
and the density of traps is small, every transfer is statisti- 
cally similar. The fast algorithm performs only one transfer 
for each cloud of electrons, then multiplies its effect by the 
number of transfers it will see. This can be implemented 
quickly in practice by sweeping one pixel's worth of traps 
up the CCD (rather than all the electrons down the CCD). 
This approach is very powerful: we confirm that it decreases 
runtime by a factor of ~ 1000 and even still allows for sec- 
ondary charge capture, whereby a trailed electron can be 
subsequently recaptured and retrailed. 

The limitation of this algorithm is that all of the cap- 
ture (and recapture) of electrons is implemented at the level 
appropriate to the size of the electron cloud in the raw im- 
age. In the real readout process, as electrons are gradually 
removed from an image peak, the electron cloud shrinks 
and fewer are subsequently captured. Similarly, as electrons 
build up in a trail, the cloud grows and becomes exposed 
to more traps. This effect is illustrated as the difference be- 
tween the dashed and solid lines in figure [T] and is most 
severe for faint sources - from which the fast algorithm even 
makes it possible to trail more electrons than are available. 
We propose a compromise between speed and accuracy by 
using one transfer to represent the first n t transfers, then 
performing a new transfer to represent the next n t and so 
on. Thus the height of local maxima slowly reduces and the 
height of trails slowly increases. We find that n t = 140 still 
provides a factor of ~70 speedup, while producing a trailed 
image within 1 electron of that produced by the full algo- 
rithm everywhere on the detector for model parameters ap- 
propriate in early 2010. 

Only one iteration of t he jBristowl (|2003al ) algorithm 
(see table 1 in lMassev et ql.ll2010r ) was required to correct 
the circa 2004 COSMOS survey. This was because the low 
density of charge traps implied only a small, perturbative 
correc tion. To correct more recent data, lAnderson fc Bedinl 
(|201Gh implemented five iterations. To qualitatively justify 
the number of iterations, it is merely necessary to test for 
convergence by calculating the difference to the corrected 
image after each step. Typical science images from early 
2010 change by only one electron in a handful of pixels after 
three iterations, and by less than an electron in every pixel 
after four. Since each iteration has a large price in run time, 
we shall henceforth stop at the third iteration. 
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Figure 2. Measured density of charge traps in the ACS/WFC detectors, as they have accumulated over the lifetime of the camera. 
Measurements assume three trap species in a ratio of 1 : 3.38 : 2.85, with characteristic release times as described in the text. Grey 
(black) points indicate survey imaging acquired with a commanded gain setting of 1 (2), and all errors are la. Separate fits are shown 
to data before and after shuttle Servicing Mission 4, plus (noisier) fits to shorter periods in grey. Hatched regions indicate times when 
ACS was offline. Points with dotted error bars show the total absolute density of traps after correction. 



3 IMAGE CORRECTION 

We use our updated CCD readout model to correct science 
imaging thro ughout the lifeti me of ACS, following the same 
procedure as lBristowl (|2003al ). The points with dotted error 
bars in figure [2] show the effective density of charge traps 
after correction, which are a factor of 20 lower than in the 
raw data and consistent with image quality in the first six 
months of operations. For the sake of clarity, equivalent post- 
correction measurements are not shown for earlier epochs, 
but th ese recover about th e same factor of 10-f5 correction 
seen in lMassev et al . (2010). Thus, ironically, as the CT1 has 
got worse, the trailing has become easier to measure and the 
correction has become more accurate! 

Figure [3] shows a region of a typical exposure, which 
was intentionally not used when measuring parameters of 
the readout model. The charge trailing that is now readily 
apparent in visual inspection of recent ACS images is suc- 
cessfully removed by our correction scheme. 



4 DISCUSSION 

We have developed a physically-motivated model of the 
readout and Charge Transfer Inefficiency (CTI) in the 
ACS /WFC detectors throughout their lifetime. We find that 



there are approximately 1.3 charge traps per pixel in mid- 
2010, split between three different species. The extended 
trails produced by these traps can be accurately modelled 
as a sum of three decaying exponentials. We also used our 
model to correct images, reducing the amount of trailing by 
a factor of ~20, to a level seen in the first six m onths of 
orbital operations. As with lChiaberge et all (2009), we still 
find no evidence for significant serial CTI (trailing perpen- 
dicular to the main trails, created by charge traps in the 
serial readout register), and therefore ignore this effect. 

When building our model, we adopted the best available 
measurements from science imaging (which we perfo rmed) 
and dark exposures (from lAnderson fc Bedinl |2010| ). The 
dark exposures were particularly useful to constrain the ex- 
tended shape of trails out to ~100 pix els and thus provid e 
better correction of object photometry dRhodes et alS^Qld ). 
Where measurements disagreed, the data's support for our 
physical model encouraged us to first extrapolate measure- 
ments obtained in the most reliable regime. 

Removing the final few percent of CTI trails might re- 
quire detailed investigation of such disagreements. In partic- 
ular, there is mounting evidence that trails behind sources 
of different flux may change in shape as well as ampli- 
tude. A slight steep ening of faint trail was also present in 
iMassev et al\ ((2010), but ascribed to uncertain background 
level and read noise. Read noise is added to an image after 
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signal induced the reado ut direction. Extrapo lating from the 
trap characterisation of iRhodes et al\ l|201(t ). we estimate 
in mid-2010 a mean shear of ~5% in galaxies detected at 
a S/N of 10. Similarly, we expect a value twice as bad for 
a galaxy at the chip gap (but zero at the edge), and about 
half as bad in a galaxy one magnitude brighter. Verifying 
this in pra ctice would require a new survey similar in size to 
COSMOS jScoville et az]|2007h . 

Most dramatically, the charge trap density increased 
~ 80% more than expected between the failure of ACS in 
January 2007 and its resumpion of activities after shuttle 
Servicing Mission 4. It is not yet clear whether this degra- 
dation is related to the decrease in operating temperature, 
the increase in temperature while ACS was offline, or coin- 
cidentally due to the ending of the solar cycle. Our current 
analysis uses almost all the suitable archival data currently 
available without yielding conclusive evidence. To resolve 
this issue, we plan continued monitoring for a further year. 




Figure 3. A typical _raw ACS/WFC science exposure from early 
2010 (HST-GO-11689, PI: Renato Dupke) before (left) and after 
(right) CTI correction. The 380 X 820 pixel area selected is furthest 
on the detector from the readout register, and the logarithmic 
colour scale is chosen intentionally to highlight the CTI trails. 



CTI, so creates spurious faint peaks that are not trailed, 
and act to spuriously steepen the true mean trail when they 
are accidentally included in the average. A physical effect 
that we do not model, but which might also affect faint 
trails, is the breakdown of the volume-d riven charge packe t 
model at very low flux levels discussed bv lShort et al\ (|201dh . 
However, while this is important in Time-Delay Integration 
(TDI) mode observations (and potentially dark exposures), 
it is not so in science imaging where a large zodiacal sky 
background is always present. If anything, the effect would 
also predict shallower trails behind faint sources, from high- 
t traps. A second physical mecahnism by which the trail 
could change shape could be the onset of surface full well 
traps above a certain flux. However, this explanation seems 
unlikely at a value of 20,000 electrons (c./. > 80,000 full 
well depth), and because such traps would have been present 
since manufacture, while almost all appear to have accumu- 
lated over time at the same rate. We shall therefore continue 
to use a single trail profile, but recommend further testing 
of this apparent shape change, for example in combination 
with mean-variance measurements at a range of flux levels 
to determine whether the shape change is gradual or discreet 
and, if discreet, whether it coincides with other discontinu- 
ities. 

The 4C decrease in the operating temperature of the 
ACS/WFC detectors in July 2006 did not affect the density 
of charge traps, or the amount of flux lost from a source. 
Howver, it lengthened their release times and the amount 
of spurious flux in the first 9 pixels behind a source fell 
by 22%, which benefits some astronomical measurements. 
Weak lensing measurements suffer by way of a spurious shear 
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